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1. Introduction 

In this work, we compute the perturbative value of large Wilson loops up to two loops, using 
the clover (SW) and overlap fermions. Using the perturbative values of Wilson loops of infinite 
length, we evaluate the shift of the b-quark mass. The perturbative values of Wilson loops of 
infinite time extent lead us to the evaluation of the static potential. For the case of clover fermions, 
we also compare our results with established results. 

The calculation of Wilson loops in lattice perturbation theory is useful in a number of ways: 
a) It leads to the prediction of a strong coupling constant fl^(mz) from low energy hadronic phe- 
nomenology by means of non-perturbative lattice simulations (for a list of relevant references, see 
our longer write-up [[I]]), b) It is employed in the context of mean field improvement programmes 
of the lattice action and operators, c) In the limit of infinite time separation, T oo^ Wilson 
loops give access to the perturbative quark-antiquark potential. Furthermore in the limit of large 
distances, /? ^ oo, the self energy of static sources can be obtained from the potential, enabling 
the calculation of (m^) from non-perturbative simulations of heavy-light mesons in the static 
Umit [g]. 

2. Calculation of Wilson loops 

The Wilson loop W, around a closed curve C, is the expectation value of the path ordered 
product of the gauge links along C. W {R, T) denotes a rectangular Wilson loop with dimensions 
RxT. There are two Feynman diagrams involving fermions, contributing to W (/?, T) at two loops, 
as shown in Fig. 1. 
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Fig.l: Two-loop fermionic diagrams contributing to W {R,T). 

The grid-like square stands for the 2-point vertex of W {R, T), whose expression is: 



W {R,T) 



-E 



ilTl) 



■Al ik)Ai{k') 5"''5{k + k') 



25^,vS (^^,/?)£sin2 {kpaT/2) 



+ 25^,vS (A:^,r)£sin2 {kpCiR/l) -4S {k^,R)S{ky,T)^m {k^a/l) ^m{kya/l) 



(2.1) 



where: S{k^,R) = sin^ [Rk^a/l) /sin^ {k^a/2). 

The involved algebra of the lattice perturbation theory was canied out using our computer 
package in Mathematica. The value of each diagram is computed numerically for a sequence of 
finite lattice sizes. Their values have been summed, and then extrapolated to infinite lattice size. 
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3. Calculation with Clover Fermions 

The fermionic part of the action contains an additional (clover) term, parameterized by a co- 
efficient, csw, which is a free parameter in the present work; csw is normally tuned in a way as to 
minimize ff{a) effects. The perturbative expansion of the Wilson loop is given by the expression: 



W{C) /N=\-Ww 



Wnlo- ^ ^'' NfXU^ 



(3.1) 



where: Wlo and Wmlo are the pure gauge contributions with the Wilson gauge action, which can 
be found in Ref. and: 



X —Xw+ X'^csw + ^Iw^iw • 



(3.2) 



We have computed the values of Xw, X^^y, and We compare our results with those of 
Ref. [|] (only Xw)\ and Ref. [|] (for m = 0). Our results for X^, X'^^, and X^^ as a function of 
mass for square loops, are shown in Figs. 2-4. 
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Fig. 2: Xw for loops Lx L. Top line: L ~ 16, 
bottom line: L = 1. 



Fig. 3: X^^^ for loops Lx L. Top line: L—1, 
bottom line: L = 16. 
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Fig.4: X^^^ for loops Lx L. Top line: L = 16, 
bottom line: L = 1. 



Fig. 5: Xoveriap for loops Lx L. Top line: 
L — 1, bottom line: L = 16. 



^In order to compare with Ref. 
errors stemming from that data. 



, we deduce their values of from the data presented there and we estimate the 
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4. Calculation with Overlap Fermions 

The fermionic action now reads 



'5f=LLv^/£»n(x,3')V^/. (4.1) 

with: Dn = Mo [l +X{X^X)-^I'^] , and: X = D^-Mo\ the index / runs over Nf flavours. Here, 
Dw is the massless Wilson-Dirac operator with r = 1, and Mq is a free parameter whose value must 
be in the range < Mq < 2, in order to guarantee the correct pole structure of Dn. 

Details on the propagator and vertices of Sp can be found in our Ref. [Q] and references therein. 
The perturbative expansion of the Wilson loop is given by the expression: 

W{C)/N=l- Wwg^ - (wnw - ~^^ NfXo,erlap^ g\ (4.2) 

where Xovertap are the values we compute. Fig. 5 shows Xoveriap as a function of Mq for square 
Wilson loops. Tables of our numerical values are presented in our Ref. 

5. Calculation of the b- quark mass shift 

In perturbation theory, the expectation value of large Wilson loops decreases exponentially 
with the perimeter of the loops: 

{W{R,T)) r^exp{-2dm{R + T)). (5.1) 
Following Ref. |^, the perturbative expansion for {W {R,T)) is: 

{W{R,T)) = l-g^W2 (/?, T) - / W4 {R, T) + ff (/) . (5.2) 
Using the expectation value of W {R, T), we obtain the perturbative expansion for 5m: 
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2(R + T) 



g' W2 {RJ)+g'[ {R, T) + {R, T 



(5.3) 



W2iR,T) involves only gluons and: W4{R,T) = {R,T) +w{ {R,T), where Wf (RJ) is the 
contribution in the pure gauge theory and {R, T) is the fermionic contribution. 

To evaluate the effect of fermions on the mass shift, we must examine their contribution in 
the limit as 7?, T — > oo. To this end, we note that our expression assumes the generic form (modulo 
terms which will not contribute in this limit): 



/ -^sin2(;,vr/2)sin2 [p^R/2) 
J [2%\ 
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Gip), (5.4) 



where p^ = 4^sin^ {Pp/'^) and G{p) ~ p^. As R,T ^ oo, the above expression becomes: 
P 
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(2 



where: p = {pi,P2,P3,(^) (for /X = 4 or v = 4). 
The fermionic contribution takes the form: 



(5-5) 



W{ {R, T) = {N^ - 1) Nf {R + T) V, (5.6) 

where V = V\y + V^^csw + VsW^lw clover fermions, and V = VoveHap for overlap fermions. The 
values of Vw, V^^^, V^^r and Voveriap have been calculated in the present work. Figs. 6 and 7 show 
their values as a function of bai*e mass for clover fermions and of Mq for overlap fermions: 
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Fig. 6: Vw, V^^, and V^^, as a function of m. 



Fig.7: Voveriap as a function of Mq- 



At one-loop order the b-quark mass shift (for N = 3) is given by: 

a5m~0.16849g^ + ^(/), (5.7) 

Using our results arrive at the two-loop expression for 5m. We list below some examples: 
The general form of 5m is {oCq = /An): 



a5m~2.1173ao + 



11.152- — — '-^V 



IN 



(5.8) 



For clover fermions, setting m = 0.0, we find: 



y = 0.00134096(5) -0.0001641(l)c5w+0.00185871(2)c|^. 



(5.9) 



These numbers agree with the result given in Ref. [@], for the case csw = 0, within the precision 
presented there. For overlap fermions, setting Mq = 1-4, our result is: 



Vov.w«p =0.000571(1). 



(5.10) 
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6. Calculation of the Perturbative Static Potential 

The static potential is given by the expression: 

aV{Ra) = -Y,ra^^^^^l^ = V,{R)g' + V2{R)g' + .... (6.1) 

where: V\ {R) is a pure gluonic contribution V2 (R) contains a gluonic part Vg (R) which can be 
found in Ref . ^ and a fermionic part F {R) : 



ViiR) 



N 



-NfF (R) . 



(6.2) 



We compute F (R) for clover and overlap fermions. For clover: 



F {R) = Fciover (R) = Fw (R) + (R) Csw + F^^ (R) 4^ , (6.3) 

and for overlap: F (R) = Foverhip {R)- The values of % F/^^ (^)' ^sw (^) Foveriap (R) for 
specific mass values, can be found in the Tables of our Ref. [||]. We also present them in Fig. 8 (for 
Csw = 1-3) and Fig. 9, as a function of R. 
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Fig. 8: Fciover as a function of R (csw = 1-3). 



Fig. 9: Foveriap as a function of R. 
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